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Abstract 



There has been renewed interest in the exploitation of Barta's configura- 
tion space theorem (BCST, fl937)) which bounds the ground state energy, 
< Egr < Supx^^^^^, by using any ^> lying within the space 
of positive, bounded, and sufficiently smooth functions, C. Mouchet's (2005) 
BCST analysis is based on gradient optimization (GO). However, it over- 
looks significant difficulties: (i) appearance of multi-extrema; (ii) inefficiency 
of GO for stiff (singular perturbation/strong coupling) problems; (iii) the 
nonexistence of a systematic procedure for arbitrarily improving the bounds 
within C. These deficiencies can be corrected by transforming BCST into 
a moments' representation equivalent, and exploiting a generalization of the 
Eigenvalue Moment Method (EMM), within the context of the well known 
Generalized Eigenvalue Problem (GEP), as developed here. EMM is an alter- 
native eigenenergy bounding, variational procedure, overlooked by Mouchet, 
which also exploits the positivity of the desired physical solution. Further- 
more, it is applicable to hermitian and non-hermitian systems with complex- 
number quantization parameters (Handy and Bessis (1985), Handy et al 
(1988), Handy (2001), and Handy, Msezane, and Yan (2002)). Our analy- 
sis exploits various quasi-convexity/ concavity theorems common to the GEP 
representation. We outline the general theory, and present some illustrative 
examples. 
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I. INTRODUCTION 



A. Deficiencies of Barta's Configuration Space Formulation 



The recent work by Mouchet (2005) develops gradient optimization strategies for im- 
plementing Barta's (1937) eigenenergy bounding procedure for the (bosonic) ground state 
energy, Egr. In particular, for any stricly positive trial function, \E' > 0, within the set 
of positive, twice differentiable, and bounded functions, C, the infimum and supremum, as 
defined below, generate lower and upper bounds to the ground state energy: 



The restriction to the ground state is because of the well known theorem that in configuration 
space, the ground state wavefunction is positive, > 0. 

Mouchet 's analysis assumes that one is working with closed form expressions for the trial 
wavefunction. Although this approach is quite fiexible, it is an incomplete resolution of an 
important, well known, deficiency of Barta's formalism (i.e. "Barta's deficiency"). Specif- 
ically, given bounds generated from an initial trial function, how does one systematically 
improve upon this, to arbitrary tightness of the bounds, within the set C? This is particularly 
important since the selection of suitable trial functions is not an intuitive process. To begin 
to answer this, within Mouchet's formalism, one requires infinite parameter representations 
for all of C. This is not possible in configuration space. 

Besides this limitation, Barta's configuration space formulation also presents significant 
multi-extrema complications in sampling the ratio ^^^j^ over all x values. Clearly, any 
alternate approach that convexifies this problem (or its equivalent), thereby requiring the 
determination of one global extremum, would be a tremendous improvement. There are 
many important problems in physics where the appearance of multi-minima, and strategies 
for circumventing these, is a major concern. One important class of methods for doing this 
simulated annealing (Kirkpatrick, Gelatt, and Vecchi (1983)). Related methods such as 
iterative annealing have impacted more contemporary research such as those studying the 
protein folding problem (Thirumalai and Hyeon (2005)). 

An additional difficulty with Mouchet's gradient approach is that it may not be the 
ideal strategy for dealing with stiff systems (such as those associated with singular pertur- 
bation/strong coupling interactions) for which very small integration steps may slow the 
global search for the infimum and supremum. (By way of contrast, the approach presented 
here does not require a gradient search. Instead a linear programming based bisection ap- 
proach proves highly effective.) 

B. A Remedy to Barta's Deficiencies: The Eigenvalue Moment Method 

In the 1980's these deficiencies were well known to D. Bcssis and his group at Saclay, 
particularly as a consequence of Barnsley's (1978) earlier studies on Barta's theorem. Their 
objective was to develop tight (converging) bounds for the notoriously difficult quadratic 
Zeeman effect for hydrogenic atoms in superstrong magnetic fields (QZE). In particular. 
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Le Guillou and Zinn- Justin (1983) emphasized an intricate order dependent, hypcrvirial, 
conformal analysis which gave numerical predicitions for the QZE ground state binding 
energy. The QZE is an example of a singular perturbation/strong coupling problem (Bender 
and Orzag (1978)), for which different methods (i.e. variational, numerical, analytical, 
etc.) can yield widely varying results. Bessis was interested in developing positivity based 
alternatives to Barta's theorem by which to asses the accuracy of of Le Guillou and Zinn- 
Justin's results. 

Independent of these concerns. Handy (1984) discovered that certain well known theo- 
rems within the classic mathematical literature known, collectively, as the Moment Problem 
(Shohat and Tamarkin (1963)) could be used to quantize physical systems through the 
generation of (geometrically) converging, lower and upper bounds, to the ground state en- 
ergy. The Moment Problem is concerned with the necessary and sufficient conditions the 
power moments of a function, Hp = J dx x^'^^x), must satisfy in order to conclude that 
the underlying function is positive, ^ > 0. More generally, Handy's interest in moments' 
quantization formulations originated from earlier studies that suggested their effectiveness in 
studying the multiscalc dynamics of certain singular perturbation-strong coupling problems 
(Handy (1981)). That is, a moments representation implicitly defined a multiscale hierar- 
chy of sensitivity to smaller and smaller scale structures. This was a precursor to wavelet 
theory (Grossmann and Morlet (1984) and Daubechies (1988), and was used, more recently, 
to incorporate continuous wavelet transform theory into quantum mechanics (Handy and 
Murenzi (1997,1998,1999)). 

Handy's eigenenergy bounding procedure exploited the confiuence of several, hitherto, 
separate results. These were: (i) the (multidimensional) Schrodinger equation with rational 
fraction potential could be readily transformed into a moment equation recursion relation 
for the power moments of the wavefunction, iJ>p — J dx xP'${x), involving the energy as 
a parameter, E; (ii) the (multidimensional) bosonic ground state wavefunction must be 
positive, "^gr > 0; (iii) the moment problem positivity theorems define an infinite hierarchy 
of constraints on the power moments, and in turn, on E. The particular positivity theorem 
used by Handy (1984) was the well known nesting constraints for the (diagonal and off- 
diagonal) Pade approximants (generated from the //p's) of the associated Stieltjes integral 
for (Baker (1975)). 

In their first collaboration, Handy and Bessis (1985) recognized that the Stieltjes-Pade 
theorems would not be extendable to multidimensions. They proposed to exploit the moment 
problem positivity theorems based on the nonlinear, Hankel-Hadamard (HH) determinantal 
inequalities. These could be extended, in principle, to multidimensions; however, this proved 
too costly. In a subsequent breakthrough (Handy et al (1988)), they realized that the HH 
nonlinear formulation could be transformed into an equivalent linearized version, which was 
then amenable to linear programming analysis (Chvatal (1983)). Handy devised an efficient 
(bi-section) method referred to as the Cutting Algorithm, which led to the generation of tight, 
converging, bounds to the QZE problem, confirming the results of Le Guillou and Zinn- Justin 
(Handy, Bessis, Sigismondi, and Morley (1988)). The entire procedure is referred to as the 
Eigenvalue Moment Method (EMM). 

Thus, historically, the EMM approach was specifically invented to bypass all of the 
unattractive features of Barta's configuration space theorem. Mouchet overlooks this in his 
review of variational (bounding) methods, despite the fact that the EMM procedure is an 
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affine map invariant, variational procedure (Handy and Murenzi (1998)). That is, the EMM 
bounds automatically sample over all affine map transformations (translations, scalings, 
rotations, etc.) of the trial functions (i.e. polynomials). Affine maps are at the heart 
of fractals (Barnsley (1988)) and wavelet transform theory (Grossman and Morlet (1984), 
Daubechies (1988)). Fractals and wavelets represent important representations for deahng 
with systems with significant multiscale structures. For this same reason, the EMM bounds 
are very good for dealing with (stiff) singular perturbation type systems with signficant 
multiscale dynamics. 

C. Moment Problem Reformulation of Barta's Theorem: The Generalized 

Eigenvalue Problem 

Despite the successes of the EMM procedure, Mouchet's work has inspired the author 
to develop a Moment Problem counterpart to Barta's configuration space theorem. Sev- 
eral new results ensue from this. The first is that the new formulation does not require 
the introduction of a moment equation. Instead, we are able to generalize the underlying 
EMM philosophy in a manner more in keeping with operator (matrix) theory. Whereas the 
moments in the EMM formulation are constrained by the moment equation relations, in 
the new formulation, all of the moments are unconstrained relative to one another. The 
energy parameter does not explicitly appear. Instead, one studies the extremal eigenval- 
ues of the corresponding Generalized Eigenvalue Problem (GEP) (Boyd and Vandenberghe 
(2004)) defined according to HT^ = AUT^, where H represents the moment representation 
operator (matrix) for the Hamiltonian, and U represents a positive definite matrix operator 
(the Hankcl moment matrix) . We arc not interested in the generalized eigenvectors , V, but 
make reference to them for clarity. Instead, it is the extremal eigenvalues, Xmin < A < Xmax, 
that are of interest, since these will directly relate to the infimum and supremum expressions 
in Barta's configuration space theorem. 

We will investigate two different versions of the above GEP system. In the first case, we 
will assume that the moments, Hp, used in defining the (H, U) operator pair correspond to 
a positive trial function (i.e. either we have a closed form for the function, \E', and its power 
moments, /ip] or we only know the /^p's, but have no closed form expression for "if). We 
can then generate monotonically converging sequences to both the infimum and supremum: 

^"l^fxi^^^ < ■ ■ ■ < Xjnin;n < ■ ■ ■ < A^m;! ^-^d Xmax;l < ■ ■ ■ < Xjnax;n < ■ ■ ■ < Supx(^^^- In 

principle, this already defines one clear advantage, since any multi-extrema features of the 
ratio ^ are circumvented by our ability to generate monotonically converging sequences. 

The (H, U) matrix operator pair, although of infinite dimension, can be studied in terms 
of their finite dimensional, upper left hand, submatrices. For each of these submatrices, we 
can define a finite dimensional, convex moment variable space, Un- In our second study of 
the GEP formulation, we will determine the optimal values, over Un, for each of the extremal 
GEP eigenvalues (i.e. Sup^^^Uni^min-Alj)) and Inf^,(,u„{,^max;n{lj)))- These wiU generate con- 
verging bounds to the ground state energy: Infi^eUniXmax-AlJ')) < Egr < Sup^^Un{^mvn;n{lj))- 
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D. Technical Preliminaries 



In order to clarify the notation used in connection with the new contributions of this 
work, we develop them in the context of a short technical overview of the Eigenvalue Moment 
Method (EMM). The EMM will also play an important role in proving some of the theorems 
introduced in this work; therefore, its review, here, will facilitate the overall understanding 
of our new results. 

1. The Moment Equation 

All of our results are predicated on being able to transform the Schrodinger equation into 
a Moment Equation representation. This is always possible for multidimensional systems 
with rational fraction potentials. For problems not of this type, it still may be possible to 
identify coordinate systems in which the transformed Schrodinger equation involves function 
coefficients that are of rational fraction form. This was the case for the quadratic Zeeman 
effect, as previously cited, for which parabolic coordinates led to the identification of a mo- 
ment equation. Despite these limitations, these restrictions still define a large and important 
class of physics problems. 

Because our intent is to develop the underlying theory, we have chosen to limit all 
discussions, and examples, to one dimensional systems, for simplicity. 

For one dimensional systems, the moment equation corresponds to a recursive, linear, 
homogeneous, finite difference equation of order 1 + m^, wherein all of the moments are 
linearly dependent on the first 1 + moments {/Iq, . . . ,f^ms}- The latter will be referred 
to as the initialization moments, or the missing moments. Once any suitable normalization 
prescription is adopted, the moments can be written in terms of the unconstrained initial- 
ization moments. The form of the normalization prescription will vary, depending on the 
nature of the problem (i.e. Stieltjes, Hamburger, etc.). In the Hamburger case, the even 
moments must be positive, the odd moments can have arbitrary signature. We want to 
choose a normalization prescription that automatically bounds the initialization moments. 
For the Hamburger case, one must have = even, and we can take fiQ + fi„is = 1- It then 
follows that \ne\ < /xo + fJ^ms = 1, for < £ < m^. Since /xq = 1 — /^ms) of the remaining 
initialization moments are unconstrained. We can now write the moment equation as 

lip = Me(p, 0) + E Me(p, P > 0, (2) 

where the Me{Pi ^) are known coefficients, nonlinearly dependent on the energy, E. We note 
that the (unconstrained) initialization moments must lie within the nis dimensional cube: 
(/xi,...,/.^J e (-l,l)-^ 

2. Stieltjes- Pade Positivity Quantization 

For parity invariant systems, the change of variables x — y^, y >0, transforms the Ham- 
burger moment problem (involving functions on the entire real axis, 3f?, fip = /J^^ dx x*'^(x)) 
into a Stieltjes moment problem (functions restricted to the nonnegative real 
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/q°° dy yP^{y)). For some of these systems, such as the harmonic oscillator problem (i.e. 
— + = E'^i^x) ), the moment equation's order becomes imity, or m^, = 0. In such 

cases, the even order Hamburger moments of the wavefunction become Stieltjes moments 
of the modified wavefunction $(|/) = — : A*2p = Up = dy y''^{y). Because = 0, 
the Stieltjes moments for the ground state become (known) nonlinear functions of E, the 
energy variable ( i.e. itp+i = Eup + 2p(2p — Uq = 1, Ui — E, U2 — 2 + E'^, etc.). 

Since the Stieltjes moments for the $-ground state are known functions of E, these then 
also determine the Fade approximants (Baker (1975)), [M\N]{^e;s)-i for the associated Stieltjes 
integral, I{s) = dy 

It is a well known theorem that if a Stieltjes measure is positive, then the [M|M] and 
[M — l\M] Fade approximants must satisfy a nested structure: 

[M-l\M]^E,;s) < [M\M+l]^E,;s)<I{s) < [M + 1\M + l]^E,;s) < [M\M]^E,;s) (3) 

Handy (1984) discovered that one could use this nested behavior to quantize the ground 
state energy, through converging lower and upper bounds. Thus, for arbitrary E, one 

generates the first Q Stieltjes moments (and all the Fade approximants that can be generated 
from them), and determines the energy interval, {Eq, Eq), of feasible energy values that lead 
to Fade approximants satisfying the above nested structure. The endpoints of the feasibility 
energy interval become the numerically generated lower and upper bounds to the ground 
state energy: Eq < Eg^ < Eq. The entire process is repeated at the next higher order 
(Q Q + 1), resulting in a reduction of the feasibility energy interval. In this manner, 
geometrically converging, lower and upper bounds are obtained. 



3. Hankel-Hadamard Determinant Positivity Quantization 



Since Fade approximants could not be easily extended to multidimensions, an alter- 
nate equivalent to the above Moment Froblem quantization procedure was required. The 
standard Moment Froblem positivity constraints, for a nonnegative function, f{x) > (ex- 
cluding distribution type expressions with zero measure support), are generally expressed in 
terms of the Hankel-Hadamard (HH) determinantal, inequality constraints, given in Eq.(6). 
These are derived through the quadratic form integral expression 



dx(j2Cr.x'^) f{x)>0, 
n=0 



(4) 



or 



/ 1^0, 1^1, ■ ■ ■ , I^N \ 
V/^AT, flN+l, ■ ■ ■ , A*2Ar/ 



|C) > 0, VC 7^ . 



(5) 



The real and symmetric Hankel moment matrix (defined in Eq.(5)) is therefore positive 
definite, with positive eigenvalues. Thus all of its subdeterminants of the following type, 
must be positive: 
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A,n,n(/i) > 0, m^even) > 0, 71 > 0, (6) 

where 



Am,n(A*) 



Det 



(7) 



For the Hamburger Moment Problem (i.e. the moment constraints leading to a positive 
function on 3?), it is sufficient to to require that A^,i(^/i^ > for m = 0, n > (Baker 
(1975)). 

It is clear that the HH inequalities are necessary conditions for any positive func- 
tion. That they arc sufficient for establishing the positivity of the underlying func- 
tion can be motivated as follows. One can approximate the gaussian, dirac distribu- 

tion, in terms of the quadratic form expansion. That is ^g^e = 2,3 j ~ 

■:(j2'i=o^(— ^^o^^ ■ Thus, for sufficiently large J and small (3 values, the HH con- 



f3V^\'^3=0 j\\ 213 

straints are essentially sampling the local behavior of a (bounded, asymptotically decaying) 
function, and requiring that it be positive. 

The HH determinants can be extended to multidimensions, as developed in the work by 
Devinatz (1957). 

The first use of the HH, Moment Problem (MP), positivity theorems to quantize the 
bosonic ground state energy was published by Handy and Bessis (1985). We briefly outline 
the essentials of this work. Through the Moment Equation's structure, Eq.(2), the moments 
are explicitly dependent on the energy variable, E, and the (unconstrained) initialization 
moments, {/xi, . . . , Hms)- So too arc the HH determinants, Ao,„(/x) = Ao,n(A*i, ■ ■ ■ , A^mJ E), 
for n < 00. Given an arbitrary E value, and (even number) moment expansion order, 
Q < 00 (thus all the moments fJ^pKQ are generated) one determines if there exists an m^- 
dimensional,initialization moment solution set, Uq-e C (—1,1)™', satisfying all the corre- 
sponding HH inequalities (Ao,o(a*) > 0, Ao,i(/x) > 0, . . . , Ag q (/x) > 0). If this solution set 
exists {Uq-e 7^ 0), then the associated E value is a possible physical ground state value, 
to order Q. If not {Uq-e — 0), then the chosen E value is not a possible physical value 
for the ground state energy. In this manner, a feasibility energy interval is determined, 
as before, in the Pade case. It can be shown that Uq e must be a convex set, if it exists 
{Ee{E^,E^) « UQ.,Ey^0). 



4- Variational- Linear Programming, Moment Problem Quantization 

Although the previous HH-MP procedure created greater flexibility in extending the un- 
derlying positivity quantization philosophy to more systems, its structure (i.e. the nonlinear 
dependence on the moments) made it too difficult for multidimensional systems. 

Any bounded, convex, set (with nonUnear boundaries), such as the Uq-e^s can be rep- 
resented as the intersection of infinitely many bounded polytopes (convex sets with hyper- 
planes as boundaries). One realizes that Eq.(5) defines this equivalent, alternative linear 
representation for the HH inequalities. That is, instead of working with a finite number of 



7 



(large dimensioned) nonlinear inequality relations (i.e. the HH constraints), one can work 
with an equivalent set of infinitely many, linear, constraints. This provided the theoretical 
breakthrough in facilitating the implementation of Moment Problem positivity quantization 
strategy. 

In order to capitalize on this linearized equivalent, one must devize an optimization 
strategy to essentially determine the optimal C's, or cutting vectors (refer to Eq.(5)). This 
required a clever combination of the moment equation formalism with linear programming 
theory (Chvatal (1983)). Thus, through the ensuing Cutting Algorithm, devized by Handy, 
given an arbitary E value, one rapidly "cuts-up" the starting (normalization) polytope ( 
the hypercube, (—1, 1)"*^) into either the null set (thereby estabhshing that Uq-e — 
and E is unphysical), or into a polytope, V D Uq-e, containing an initialization point, 
jl =jjli, . . . ,jlm^) for which all the associated Hankel matrices are positive (concluding 
that pL G Uq-^Ei hence Uq-^e 7^ 0, and that particular E is, & possible value for the ground 
state energy). 

The above, entire, procedure is the Eigenvalue Moment Method (EMM), which was used 
to solve the previously cited quadratic Zeeman problem. 
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II. MOMENT PROBLEM REFORMULATION OF BARTA'S THEOREM 
AS A GENERALIZED EIGENVALUE PROBLEM 



In this work we directly transform Barta's configuration space theorem into a Moment 
Problem representation. By so doing, we obtain a theoretically complete formalism that 
addresses and solves Barta's configuration space deficiencies. Our analysis is limited to 

Hamiltonian systems which can be represented in terms of a moment representation. As 
previously noted, this corresponds to a large, and important, class of physically interesting 
systems. 

In contrast to the EMM approach which restricts itself to the set of moments satisfying 
the physical moment equation, the new formalism makes no such restriction. Previously, 
Uq referred to a subset within the domain of initialization moments (i.e. the subset of 
initialization moment values whose Moment Equation generated moments, up to moment 
order Q, satisfy the HH positivity constraints): 

Definition (Sec. I) 



Uq = {(/xi, . . . = Me{p,0) + J:T=i ME{p,i)fie, < P < Q, and > 0, < 

n<f,}c(-l,ir^ 



Note that the elements of Uq implicitly must satisfy some, physically motivated, nor- 
malization prescription. 

The new definition of this same notation will be 

Definition (Modified) 

= {(a*o, • • . , iJ,Q)\where Ao,n(Ai) > 0, < n < f }. 

That is, Uq refers to the domain of Hamburger moments, up to moment order, Q, 
satisfying all the corresponding HH positivity constraints (and not constrained to satisfy 
any moment equation). In either case, Q is implicitly an even number. 

It is also implicitly assumed that the elements of Uq must satisfy some physically moti- 
vated normalization prescription. 

Within each of the finite dimensional subsets, Uq, Barta's relations manifest themselves 
in terms of a Generalized Eigenvalue Problem (GEP) (Watkins (2002)) 



where H and U correspond to finite, real and symmetric matrices, to be defined in Sec. HI. 

The matrix elements will be linear in the moments. Because of the restriction to Uq the 
U-Hankel matrix is positive definite. The matrices (H,U) are designated as a symmetric 
pair. As will be seen by the explicit example discussed below, whereas all of the moment 
variables in Uq contribute to the structure of H, a reduced number of these contribute to 
U; however, this reduced number still guarantee the positive definiteness of U. 

Through a Cholesky decomposition (Watkins (2002)), U = R*R, the GEP is transformed 
into a standard, symmetric matrix, eigenvalue problem. 



The extremal eigenvalues \nin < A < Xmax, satisfying the above, or alternatively 
Det{ii — AU) = 0, are defined by: 



(8) 




(9) 
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Xrmnin) = luf-^ Slfi^^lE! ^ {/^} ^ Uq. (10) 

^ (C|U(/i)|C) 

and 

XmaM = {/^} e Wq. (11) 

^ (C|U(/x)|C) 

Of course, we also have XminilA — —Sup-^ ^"^"'"'■^■'J^^ . 

C {C|u(m)|C>^ ^ 

(CjH(M)^C ) 

The above ratios can also be written in terms of ' . Let Ah*u'"''^(/^) denote 

(C\u{p)\C ) 

{C\C) 

the extremal eigenvalues of the H and U matrices, respectively. Since \ij'^{ix) > 0, if we 
also assume that A^5"'(/x) > (for simphcity), then i^e^ < \nin{^^) < XmaxifJ') < x^^^"(py 

However, if Ah*"(/x) < 0, then ^^in^^^ < X^inilA- These expressions will prove important 
later on. 

Physicists refer to any function which, locally, lies below the tangent plane as a convex 
function. Mathematicians define a function, /(x), as convex if the set {{x,y)\y > f{x)} 
is convex. Thus, what a physicist would regard as a concaved function, is referred to as 
a convex function by mathematicians. Throughout this work we will use the physicist's 
definition, except in a few cases where we cite the actual theorems (placing quotation marks 
around them), as they appear in the literature. 

Definition 

/ Concaved\ _ / "Convex" \ ^ 

\ Convex ) ^^^[p^^^^^ctsts] = I "(7oncaf erf" j ^^^I"*"*^^"*"**"""''] 

It is a well known theorem that the smallest eigenvalue of a symmetric matrix is a convex 
function with regards to the matrix elements as variables. This led to various alternative 
algorithmic strategies (i.e. gradient methods) for implementing EMM (Handy, Giraud, and 
Bessis (1991), and Handy, Maweu, and Atterberry (1996)). 

The smallest eigenvalue of the generalized eigenvalue problem is not convex as such, but 
it shares the good fortune that it does not have any multi-maxima or saddle points. That 
is, there can be regions of relative flatness. Thus, although there can only be one global 
maximum value, the points in Uq corresponding to the global maximum may not be unique. 
These properties are what the mathematicians refer to as quasi - "concaved". 

Despite this, in the infinite limit (all of moment space) the physical problem strongly 
suggests that the point, in the moment variables space, corresponding to a global maximum 
for XminilA-i unique. This is because there can only be one physical ground state (the 
ground state energy can only be associated with one point in Woo)- Since we are only 
interested in obtaining tight bounds for Eg^ these issues do not affect this objective. 

We now revert to the mathematics nomenclature which is opposite to the physicist's 
intuitive interpretation, as previously noted. 
Definition: A function, f{'jt) : W ^ is quasi- "concaved" if 

/(s7?i + (1 - s)7?2) > mm{/(7?i), /(7?2)}, < s < 1, 7?i^2 e ^- (12) 
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It therefore follows that for a quasi- "concaved" function there can be flat regions where 
the function stays constant. However, along the one dimensional path defined by < s < 1, 
there can be no local minimum. If the function is strictly "concaved", then any local 
differential search will always yield a path which leads to the global maximum. However, 
for quasi- "concaved" functions, within regions of flatness, more effort may be required to 
flnd a path that leads to a global maximum. 

We have the important theorem (mathematical nomenclature) 
Theorem # 1: XmaxilJ') is a quasi- "convex" function of the matrix elements (H,U), 
which are hnear in the moments (Siddharth (2005), Boyd and Vandenberghe (2004)) . 

Similarly, A^j„(/i) (being the negative of a quasi- "convex" function, refer to discussion 
following Eq.(ll)) is a quasi- "concaved" function. 

From the deflnition of the extremal eigenvalues one has that Amm(Ai) < ^maxi^j)- How- 
ever, there will be moment elements in IAq for which the extremal eigenvalues coincide. This 
will be the case for those satisfying the moment equation, up to moment order Q (i.e. those 
moments that also satisfy the EMM moment equation) . 

Definition: Denote by ^ = {fiQ, fii, .... fig} G Uq, an element of Uq that also satisfles the 
moment equation (Eq.(2)), for the given E value. Such a point, by deflnition, automatically 
satisfles the EMM positivity constraints. This is only possible if £■ e {Eq,Eq). That is, it 
must lie within the EMM eigenenergy bounds. Now use it to generate the (H, U) symmetric 
pair matrices, as deflned in the next section. It will be shown in the following section that 
Theorem #2 



We shall denote by Uq-emm C Uq, the subset of points satisfying the EMM conditions (i.e. 
satisfles the moment equation up to order Q, and the positivity conditions, for an E value 
that must lie within the EMM bounds). 

Let us explicitly distinguish the extremal GEP eigenvalues for each Un of dimension n + 1 
by the notation: \„iax/min:n{fJ-)- Also, let C denote the set of functions that are positive, 
bounded (exponentially decreasing), and continuously differentiable up to the second order. 
In the next section we shall show that if ^ e C, and Hp — J dx xP^{x),p < oo are its 
moments, then 
Theorem # 3: 



^minifJ'E) — E — \rnax{l^E)- 



(13) 




(14) 




(15) 



and 




(16) 



for each * e C. 
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This is an interesting result, particulary when combined with sequence acceleration meth- 
ods, since it allows one to determine Barta's lower and upper bounds for a function whose 
moments are known, even when the function is not given in closed form. Also, the monotonic 
nature of the results may prove useful in circumventing potential multi-extrema features of 
the ^ ratio, when evaluated in terms of Barta's configuration space formulation. 

Let us now define the Sup and Inf of the extremal eigenvalues over their finite dimensional 
convex domain, lAn- 

( ^7^^- ) = ( (^--;-(/^)) \ _ (17) 

From Eq.(13), since Un^EMM C we must have that A^"^.„ > , the EMM upper 
bound. Similarly, A^^.„ < E^, the EMM lower bound: 
Theorem # 4 

^IiT-f <^ ^EMM—lowerbound <^ ^ <^ ^E M M —upperbound ^ \^Sup fl8) 

and, in the infinite limit, n ^ oo: 

\Inf ^ \Inf ^ p ^ y^Sup ^ ySup /-.qn 

Important Assumption/Condition Although W„ will be a bounded convex set through 
the normalization conditions used, it is also important that its boundary {6U) not include 
points at which the positive dcfinctness is lost. That is, if ju| G we do not want 
^Qj<niMi) = 0, or Au*"(iu|) = 0, using the notation in the discussion following Eq.(ll). 
If this is were to happen, then the A^"fj.„, as well as Xm/x-ny could become singular (i.e. 
+CXD, — oo, respectively). In the present application, we can insure the above by simply 
imposing additional moment inequality constraints associated with any (rough) upper bound 
to the ground state energy. This will be clarified in the last section, where we implement 
the numerical analysis on a specific problem. 
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III. PROOF OF THEOREMS 



In this section we develop the basic relations and prove the various theorems previously 
quoted, with the exception of Theorem # 1 which is a standard result in optimization theory, 
particularly in the context of mathematical economics. We will be limiting our discussion 
to the one dimensional case, for simplicity. 

Two crucial elements are required for deriving Barta's theorem in Eq.(l). The first 
is that in the configuration space representation, the bosonic ground state wavefunction, 
'i/gr{x), must be of uniform signature, and thus can be taken to be positive, "^gri^) > 0. 
Accordingly, given any trial function, "if, of arbitrary signature, and with a bounded and 
continuous second order derivative, one obtains the zero identity for the scalar product: 
{"^grKH — Egr)\'^) = 0, where Egr is the ground state energy, and H is the Schrodinger 
equation hamiltonian. Therefore, H — Egr, when applied to must have a zero at some 
location 



The second assumption is that if the trial function is stricly positive, ^ > 0, then the 



Egr E {R{x)\\/x e 3fJ}. This leads to Eq.(l) or Inf^R{x) < Egr < Sup^R{x). 
Positivity is an important cornerstone of Barta's theorem. 

Let S and X denote the supremum and infimum, respectively, for an arbitrary trial 
function, ^f, lying within the set of functions, C, which are positive, bounded (exponentially 
decaying), and have continuous, finite, second derivatives: 




(20) 



range of the function R{x) 



must define a bounded subset of 3? that contains Ear'- 




(21) 




(22) 
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Define the configurations: 



L^{x- \i) = —-V - A; > 0, ^ A, < X, (23) 



and 

'H<if(x 



^{x; A„) = A„ - (^^) > 0, ^ A„ > 5. (24) 

Although the trial functions must be positive (i.e. strictly positive) if they are to be easily 
used within Barta's procedure, the {U^i,{x), Lq,{x)} functions can be nonnegativc (provided 
the zeroes correspond to sets of zero measure) and still generate strictly positive HH deter- 
minants. Accordingly, 

(H-Xi)^(x)>0, ^ Xi<I, (25) 

and 

(A„ - H)^f{x) >0, ^ Xu>S. (26) 
Let us focus on the first relation: 

$A,(x) ^{H- A/)*(x) > 0, A; < I. (27) 

For any \f (x) e C, ^Xi{x) must be integrable and positive almost everywhere (i.e. non- 
negative). Thus, its power moments must satisfy the standard positivity relations of the 
Moment Problem (Shohat and Tamarkin (1963)), as discussed in the Introduction. 

We are restricting our analysis to hamiltonians with rational fraction potentials, since 
these are the ones most easily transformable into a moment equation representation. For 
simplicity, the following discussion assumes that the potential is of (multidimensional) poly- 
nomial form. The generalization to singular potentials is straightforward, and briefly dis- 
cussed below. 

In order to make our analysis more transparent, we will consider the case of the quartic 
potential problem: H = —d^ +x^. Then ^x^ix) = {—dl +x'^ — Xi)'if{x), and we can generate 
the power moments of the LHS, based on those of \E'(a;). 

Define the power moments of the trial function by: Hp = dx xP^{x), p > 0. By 
assumption (i.e. ^ > 0), these must satisfy the Hankel-Hadamard (HH) determinantal 
constraints for the Hamburger moment problem: A„i^n(^l^^ > 0, for m = 0,n > 0. These 
constraints are required for all positive (more generally, nonnegative) functions on the real 
axis. 

The power moments of ^Xi{x), 

/+0O 
dx xP^Xiix), p > 0, (28) 
-oo 

satisfy (i.e. upon substituting the <l>/^ relation and performing the necessary integration 
by parts) 
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= - 1)a*p-2 + A«p+4 - Ai//p, p > 0. 



(29) 



If A; < X, then the v moments generate the Hankel matrix that must satisfy the (fifi) 
positivity constraints: 



Ao,7v(z^(A;)) = Det -{Ui + n2)(ni + - l)/ini+n2-2 + A*ni+n2+4 - A; /X, 



ni+n2 



>0, (30) 



/ 



for all < ni, n2 < A'" < oo. 

The form of the finite dimensional Hankel matrix in Eq.(30) is symbolized by H -A;U, 
with U the positive definite Hankel matrix for ^'s moments. 

We note that the {z/q, . . . , i'2n} moments used to define the Hankel moment matrix 
for depend on the {/xo, • • • ^^2N+-i} moments of \E'. Thus, we are working within the 
moment space Uq where Q = 2N + 4. Notational consistency would suggest that in the 
following discussion we make reference to \min;QM{lA where, Qn = 2A^ + 4. To streamline 
the discussion, we will simply use the notation Xmin;N- 

Define by \min;NilJ') the smallest zero satisfying (the /i dependence is not explicitly given, 
for greater clarity) 

Ao,iv(«^(Ami„;Ar)) = 0, (31) 



or 



Det(H{^) - Xrmn,NU{fi)) = 0. (32) 

Accordingly, (— oo, Ami„;Af(;u)) D (— oo,X). This is because, so long as A; < I, Eq.(30) must 
hold. Hence, any root in the A-variable domain must be larger than X. 

We now show that the sequence {\min,N{fJ^)\N > 0}, must be nonincreasing, or 
^min,N{l^) > ^min,N+i{l^) ■ The casicst way is to exploit the positive definiteness of U, 
which leads to the Cholesky decomposition U = R*R, where R is the unique, upper tri- 
angular matrix, with positive diagonal entries. Its inverse is also of upper triangular form. 
Accordingly, \min;N is also the smallest zero for the equation 

Det ( R-*HR-i - A^i„,jvl ) = 0. (33) 

Of course, this is also the smallest eigenvalue of the indicated real, symmetric matrix, which, 
in turn, define a nonincreasing sequence; thus proving the previous claim. 

With regards to the A„'s, an analogous result follows. Thus, if A^ > «S, then 

Ao,Af( - i^{K)) = Det ^ A„ Hni+n2 + {ni + n2)(ni + ^2 - l)/i ni+n2-2 Atni+n2+4 ^ > 0, 

(34) 

for all < ni, 722 < A^ < oo. 

We now define Xmax;N{lJ') as the largest root satisfying 
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^0,N{ — l^{^max;N)) — 0, 



(35) 



or 

Det{Xmax,N'i^{fJ') - H(//)) = 0. (36) 

(Note that Xmax-.N is also the largest root of the Generalized Eigenvalue Problem, Det(li{iJ,) — 
Amax,7vU(/i)^ = 0.) It then follows that {Xmax;N , +00) D (5, +00), and they form a non- 
decreasing sequence: \max;N{lJ) < Xmax;N^i{lJ)- This follows from the observation that 
N is the smallest root of 
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Det(^ — H(//) — {—Xmax,N)^{lJ')^ — 0, and through the Cholesky decomposition of the 
positive U matrix, the —Xmax-.N form a non- increasing sequence of smallest eigenvalues for 
the finite and symmetric matrix: R~*(— H)R~^. 



A. Theorem # 3 

The first part of Theorem #3 ( Eq.(14) and Eq.(15)) follow from the previous results. 
The latter part of Theorem #3 (Eq.(16)) results from the fact that the only A; values 
satisfying all of the HH positivity inequalities are those obeying A; < I, similarly for A^: 
hm„^oo Arrim;n(A*) = I, and hm 

n— >oo Am,aa;;n(/^) ^■ 



1. Alternative Derivation 



We include some additional remarks that provide a different perspective on all the above. 
We limit the discussion, for brevity, to the \min;n case. 

The HH positivity theorems require that all of the HH determinants, for the //-moments, 
be positive, as functions of A;. To simply investigate the positivity properties of one of these 
determinants is insufficient. Thus, we can either work with the quadratic form inequality in 
Eq.(5), adapted to the fp(/x, A^) moments (Eq.(29)), for the (iV + 1) x (iV + 1) dimensional 
Hankel matrix, 7iN{j^)' 

( uo{ii,Xi),ui{ii,Xi),...,un{ij>,Xi) \ 
ui{ii, Xi), U2{ii, Xi),..., lyN+iifJ', Xi) 



[C\ 



\C) > 0, VC 7^ 0, 



(37) 



V UN{fi, Xi), l^N+lifJ', Xi),..., Z/27v(Ai, Xi) J 

or we can work with the + 1 HH determinants: 

Ao,„(i/(/i,AO) >0, n = 0,...,7V. (38) 

The set of A/ values satisfying Eq.(37) must correspond to a convex set (Chvatal (1983)) 
since it represents an infinite set of linear inequalitites in the A- variable. Thus, the feasibility 
A set must be a semi- infinite interval. That is, if Eq.(37) is satisfied by the two values 

— ) 

(^1 - (ni + 712) (ni + 712- 1)a* ni+n2— 2 ~l~ A*ni+n2+4 A^ lln\+n2 1^) > 0, (39) 



(^1 - (ni + 712) (ni +712- l)Hni+n2-2 + A*ni+n2+4 



X\ ^Vni+n2|^) > 0, 



(40) 



then it must be satisfied by all A; = sX\'^^^ + (1 ~ ^)Xi'^^\ for < s < 1 (simply multiply 
each of the above two inequalities by s > and 1 — s > 0, respectively, and add); thereby 
establishing the convex nature of the set of allowed A; values. 

Since A; G (— oo,X) satisfies Eq.(37) (From Eq.(25)), it follows that there exists a 
"largest" semi-infinite interval {—00, Xi-n), with I < Xi-n, satisfying all of Eq.(37). It is 
clear that A;;jv+i < Xi-m, since from Eq.(37), the quadratic form inequalities for the TiN+iij^) 
Hankel matrix includes all of those corresponding to the T-Cn{^) case. The Xi-n value must 
then be the smallest root of the equation Det{7iN{^{lJ', Xi)) — 0. That is, Xi-n — Xmin;N- 
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2. Extension to Rational Fraction (Singular) Potentials 



The preceding, alternative, proof also shows us how to extend our results to the case of 
rational fraction type potentials. Consider the perturbed quartic potential V{x) — x^ + -^t^- 
Limiting ourselves to the "infimum" case, for simphcity, we see that Eq.(27) can be modified 
by multiplying both sides by the positive denominator polynomial, + 2: 

{x^ + 2)<^x,{x) = ( - (a;2 + 2)dl + {x^ + 2)x^ + 1 - Xi{x^ + 2))^{x) > 0, \i < I. (41) 

Thus, the R.H.S. generates a positive (Hankel) matrix, and one can procede to define the 
corresponding Xmin;N, which satisfies aU the relations described above. 

Thus, in general, as long as one multiplies Eq.(27), or its multidimensional counterpart, 
by positive "denominator" type polynomials, in order to achieve a Hankel matrix structure, 
then all of our results apply. 

3. Additional Remarks 

Given that the function ^ can have multiple extrema, having a systematic method of 
computing the infimum/supremum (as opposed to searching over all local extrema) may 
make the above results very convenient. That is, one can determine the infimum (X(/i)) 
and supremum (<S(//)) by studying the asymptotic limits of the \min;n and \max;n-i combined 
with sequence acceleration techniques, where possible. 

Another important aspect of the previous results is that we can now extend Barta's 
result to positive functions which may not be given in closed form, but whose moments may 
be known. We provide one example of this in the next section. 

B. Theorem # 2 

We now prove Theorem ^2 by way of the qiiartic potential problem. Whereas in the 
previous proofs we were working with an infinite set of numbers, {^p\p > 0}. known to be 
the moments of a positive function, we will now be working with a finite set of moments 
that satisfy the Moment Equation, as well as the corresponding positivity theorems. 

Let us assume that the {iip\0 <p < P} moments satisfy the moment equation, 

-p{p - 1)/Xp_2 + -Eiip = Q, < p < P - 4. (42) 

Alternatively, 

= E/jp + p{p - l)/ip_2, < p < P - 4. (43) 

Clearly, this recursive relation separates into the even and odd order moments. There are 
more efficient ways of dealing with such relations, however, for our immediate purposes, the 
above is satisfactory. Also, we implicitly assume that some normalization has been chosen. 

Let P = 2M. Again, we assume that the moments {hq, Hi, . . . , ij,p=2m} £ Up=2M;Emm C 
^p=2M, satisfy the moment equation and all the HH determinantal inequality conditions 
that can be generated from them, for some E value. Thus Ao,n(A*) > 0, for n < M. From 
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Eq.(30), we see that the U matrix involves all of the moments up to order //2Ar, where N is 
to be determined. The H matrix involves the highest order moment, fJ^2N+4- Thus, we want 
2N + 4 = 2M. That is, the highest dimension Generalized Eigenvalue Problem (GEP) is 
A^ + 1 = M-1. 

Prom Eq.(42), it follows that for the special set of moments being considered, we have 
H = E\J. The corresponding GEP problem becomes 

Det(ll - AU) = Det(^EV - AU) = {E - \f+^Det(v) , (44) 
where Det(\J^ = Ao,Ar(A*), revealing its (N + l)-th order degeneracy. Hence 

^min;N{fJ'E) — ^max;N{fJ'E) ■ (45) 

In summary, the GEP problem becomes extremely degenerate for those moments satisfying 
the Moment Equation, as well as all of the corresponding HH positivity constraints. The al- 
lowable E values are those generated through the EMM procedure corresponding to moment 
order P = 2M. 



C. Theorem # 4 

Define the supremum of the smallest GEP eigenvalue by Xmin;Q — S'^Pn&AQ^min;Q{lj)'^ and 
the infimum of the largest GEP eigenvalue by A^f^..^ = Inf^gu^Xmin^qilJ')- Since the EMM 
related set of moments satisfy Uq-emm C Uq, and on Uq-emm the extremal eigenvalues 
are degenerate, it follows that the EMM upper bound, must be a lower bound to A^^.g. 
Likewsize, the EMM lower bound, must be an upper bound to X^^I^.q- This confirms Eq.(18). 

In the Q — > oo limit, the entire moment space, for a given set of moments jt = 
(/io, • • • , yUj^oo), we must have that Xmin;oo{l^) < Egr < Xmax;oo{l^) , from Eq.(16). Clearly 
then A^?.^ = Egr = X^^f^.^; confirming Eq.(19). 

Again, we implicitly assumed that Uq satisfies some physically motivated normalization 
prescription. 
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IV. SOME NUMERICAL RESULTS 



We will use the quartic potential to illustrate, numerically, most of the previous results. 
The ground state energy is Eg^ = 1.060362090484. 

A. Eqs. (14-16), using *(x) = 7Ve~^^ 

This is a trivial example. One has — 2 — Ax^ + x^. The infimum is Inf{2 — 4x^ + 

x^^ — —2. The even order, gaussian function power moments, satisfy the recursion relation 

— (j^^fh-i P — 0- Normalizing according to //q = 1, determines the normalization 
factor J\f. 

Since the Gaussian function, = also satisfies —d'^'^g + Ax'^'^y = 2^^, an alternate 
recursion relation for the moments is —p{p — l)/ip-2 + 4/Xp+2 = 2/ip. Using this, one can 
transform the matrix structure in Eq.(30) (i.e. —p{p — l)Aip-2 + IJ'p+4 — ^l^p) into — 4/i,p_|_2 + 

— (A — 2) //p. The Generalized Eigenvalue Problem results, for this problem, are given 
in Table I. Note that the convergence is slow, but consistent with the various Theorems. We 
also note the curious repetitive, non-repctitivc, structure manifested by the eigenvalues. No 
sequence acceleration analysis has been attempted. 

B. Eqs.(14-16), using ^'(x) = |$(x)p, where <I>(x) satisfies the PT-invariant 
Schrodinger equation —9^$ — (ix)^$ = 

We now investigate the utility of the previous formalism when the positive trial function 
is not known, in closed form; although the moments are (numerically) known. For this 
cxccrsize, we could take ^(,t) to be the (positive) ground state of any (solvable) Schrodinger 
potential problem, —d^^s + K('i')^s('i') = E^^! g[x). In such cases, one would determine the 
Inf/Sup of (i.e. Hi = -dl + x^), ^ = -^ + x^^ E,- Vs{x) + x^. This would then be 
a trivial analysis. 

Instead, we pursue a different class of problems whose differential structure does not lead 
to an easily calculable set of Barta bounds. Such is provided by the class of non-hermitian 
systems that have received much attention in the context of PT-symmetry breaking systems 
(Bender and Boettcher (1998)). The simplest example of this is the well known V{x) = 
— {ixY system, which we write as — 9^$(x) — {ixY^{x) — £^{x). This system admits 
only real cigcncncrgics; however, its cigcnstatcs arc all complex functions. Nevertheless, the 
probability density, ^{x) = |$(x)p = S{x) > 0, satisfies a linear, fourth order differential 
equation (Handy (2001)): 

( - ^dlS{x) - A^d^S{x)) + Ax^S{x) = 0. (46) 

Although all of the bound states of this system arc positive, we shall work with the one 
corresponding to the smallest £ = 1.1562670719881133. The Hamburger moments of the 
even function, S{x), satisfy a simple recursion relation (Handy (2001)): 

4Atp+7 = (p + 4)p(p - 1) (p - 2)//p_3 + 4:Sp{p + 4)/ip_i, (47) 
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for p > 0. The GEP-moment analysis, given in Table II tells us that Inf(^^^^) < -1.57. In 
order to verify this, we can implement a Runge-Kutta analysis on S{x), in order to calculate 
— If one is not too careful (i.e. implementation of a naive second order finite differencing), 
a significantly wrong answer is obtained (i.e. Barta's infimum is 0(.3)). Instead, by using the 
relation (^^) = _£^£!±t:|^^^£:a^ + ^4 ^ 2£-2\m^ + x\ the resulting expression 



V S(x) J ^1 * 

lends itself to a more accurate Runge-Kutta verification, yielding the (approximate) Barta 
infimum as —1.782. This is very consistent with the GEP generated results in Table II. 

C. Generating Converging Bounds for Quartic Potential (Theorem # 4): 
The need for a rough upper bound to the energy 

In this last example, we will not work with a fixed set of moments for a positive 
trial configuration. Instead, we will implement an optimization procedure for determin- 
ing Sup{Xmin-^Q{iJL)^ and Inf (^Xmax-^Q{fi)y for iieUq. Contrary to the Gradient analysis in 
Mouchet's (2005) work, we can determine these quantities by combining the EMM, linear 
programming based, "Cutting- Algorithm" (Handy et al (1988)) with a bisection method in 
the A-variable space. 

1. Defining Uq 

The Uq space is defined as the set of Hamburger moments, {(/xq, . . . , A^q)}, satisfying 

01 I'^i) > 0, V'^i ^ 0, < m + 712 < Q. (48) 

In addition, for the case of Sup{\min;Q{lA) ■, we are interested in the set of A;'s satisfying 
(i.e. Eq.(30)) 

(^2! - (m + n2){ni + n2- 1)// ni+n2— 2 ~l~ A'ni+n2+4 \l^ni+n2 

1^2) >0,V^2 7^0, (49) 

< ni + n2 + 4 < g. 

Whereas, for the Inf(^Xmax;Q{l^)^, the latter set of inequalities are replaced by (i.e. 
Eq.(34)) 

02\Kl^m+n2 + + 'n2){ni +n2- l)l^ni+n2-2 " /^ni+n2+4|^2) > 0, V'^2 7^ 0, (50) 

< ni + 712 + 4 < g. 

2. Normalization Prescription: Bounding Uq 

One must also impose some normalization condition. A choice that leads to a bounded 
Uq set is: 

even) ^- 

(51) 
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To study the consequences of this, note that the physical moments for the ground state 
wavefunction (i.e. assume ^ = ^p^) satisfy 

= r dx x^'^{x) + [ dx x^^^ix). (52) 

The p — even moments must be positive and satisfy 

< IJ.p=even = dx X^*(x) + / dx X^*(x) < / dx *(x) + / dx X^'^lx), 

J-l Jx^[-l,l] J-l Jx^[-l,l] 

(53) 



or 



< Hp=even < /^O + /^Q = 1- (54) 

For the odd order moments, a similar set of relations ensues for \iip=odd\ < /-i \x^\^{x)-\- 
Sx^[-i,i\ dx \xP\^{x) < iJ,o + fJ'Q — 1- Thus, we have 

-1 < IJ'p=odd < +1- (55) 



3. Linear Programming - Bisection Algorithm for Determining A^"^ and 

The following algorithm implicitly makes use of the quasi-convex nature of \min;Q{l^) and 
the quasi-concave structure of \max;Q{lJ^) for G Uq. 

We outline the basic structure of our computational algorithm. Assume that for a trial 
positive solution, (^Uq, • • • ,/Uq), we have determined its corresponding extremal eigenvalue, 
^mm;Q(/^*)- Within the interval [Amm;Q(Ai*)) oo), we pick an arbitrary point, A^, and use 
the EMM "Cutting- Algorithm" to determine if there exists a point in IAq, satisfying the 
normalization conditions, as well as Eq.(49), for A; = Aq. There are two possibilities: 

(A) If there is such a point, then we repeat the entire procedure, but within the interval 
[A„,+oo). 

(B) If there is no such point, then the entire procedure is repeated within the interval 

['^mm;Q(A**)) '^J ■ 

The objective is to eventually generate a reducing sequence of intervals, [Xan ^a^] ^ • • • ^ 
[A(j., Aa-^J, until an acceptably small interval is attained. The endpoints will tightly bound 

\ Sup 
min;Q' 

For the \max;Qi ^ similar procedure is required. Thus, one would select a point within 
the interval {—oo, \„,,ax:QifJ'*)) ■ Upon picking an arbitrary point within this interval, Aq, 
one would then determine the existence, or non-existence of a /x-point lying within Uq, and 
satisfying the normalization conditions. Such a point must also satisfy Eq.(50), for A„ = A^. 
If there is such a //-point, then the entire procedure is repeated for the interval (— oo, Aa). 
If there is no such point, then the updated interval is (Aa, ^max;QifJ'*)) ■ 
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4- The Need for a Rough Upper Bound 



All of the above is contingent on making sure that the boundary of Uq include no 
points at which the U(/x) matrix has zero eigenvalues. This was emphasized previously. 
The adopted choice of normalization, if not supplemented with additional linear constraints 
on the moments, include such singular points. Specifically, because the U matrix only 
includes moments up to order ^2N, while the H matrix includes the additional moments 
{A*2Ar+i, • ■ ■ , A*2Af+4}, one possible boundary point could be all of the first 2N + 4 moments 
set to zero (i.e. A*o<n<2Ar+3 = 0) and the last moment set to unity, 112N+A = 1- To avoid 
these, and other such possibilities, any rough upper bound for Egj. will help in restricting 
IAq to avoid such boundaries. 

Let Epub » Egr denote a poor upper bound to the ground state energy. The true moment 
equation for the quartic problem is Egr/ip = —pip — 1)a*p-2 + A*p+4- Taking p = even, we 
have 

Epub l^p + P{P - l)/^p-2 > P = even. (56) 

Thus, these additional inequahty relations will lead to proper Uq sets. In Table III, we take 

The results in Table III confirm the all of the above theoretical results. Note that 
the EMM bounds will be, generally, tighter than those derived from a "moment problem 
extension of Barta's theorem". The calculations were done using the Stieltjes form for the 
moments (i.e. all odd order Hamburger moments were set to zero, abinitio, fj^odd = 0). The 
results in Table III confirm that knowledge of a rough upper bound for the ground state 
energy lead to converging bounds for the ground state energy. This result is similar to that 
developed in a Euclidean time reformulation of the EMM philosophy, as applied to positive 
matrices (Handy and Ndow (1992)). 



23 



TABLES 



TABLE I. Quartic potential results using Me ^ trial function; Barta's lower bound is -2 (Note 
that Dim = A + 1 and Max. Moment Order, Q, satisfy (5 = 4 + 2A) 
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-1.67629 


60 (122) 


-1.841515 


93 (188) 


-1.89958 


28 ( 58*) 


-1.68012 


61 (124*) 


-1.841519 


94 (190) 


-1.90236 


29 ( 60) 


-1.68637 


62 (126) 


-1.85045 


95 (192) 


-1.90236 


30 ( 62) 


-1.68637 


63 (128) 


-1.85045 


96 (194) 


-1.90422 


31 ( 64) 


-1.72107 


64 (130) 


-1.85818 


97 (196) 


-1.90422 


32 ( 66) 


-1.72107 


65 (132) 


-1.85818 


98 (198) 


-1.90529 


33 ( 68) 


-1.74204 


66 (134) 


-1.86361 


99 (200) 


-1.90529 










100 (202) 


-1.90576 










101 (204) 


-1.90576 



* Except for these entries, all others appear in pairs, to 20 significant figures 
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TABLE II. Quartic potential results using as trial configuration the solution 
corresponding to PT-invariant, non-hermitian, system —d'^^{x) — {ix)'^^{x) = £^{x), for 
£ = 1.1562670719881133. Barta's lower bound is approximately (-1.782) , based upon Runge-Kutta 
integration. 



Q 




Q 




4 


0.7651830316 


30 


-1.412946343 


6 


-0.3701497316 


32 


-1.412946343 


8 


-0.5797495842 


34 


-1.466405630 


10 


-0.5797495842 


36 


-1.466405630 


12 


-0.9175699949 


38 


-1.471430659 


14 


-1.025936484 


40 


-1.505086215 


16 


-1 025936484 


42 


-1 505086215 


18 


-1.202683806 


44 


-1.535124305 


20 


-1.202683806 


46 


-1.535124305 


22 


-1.202840090 


48 


-1.530286871 


24 


-1.349131584 


50 


-1.556428376 


26 


-1.349131584 


52 


-1.556428376 


28 


-1.369576335 


54 


-1.579966618 






56 


-1.5799666184 






58 


-1.5799666184 






60 


-1.5882508326 
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TABLE III. Results of the Quasi-Convexity/Concavity Analysis (i.e. "Moment Prob- 
lem Reformulation of Barta's Theorem") Applied to the Quartic Potential Problem: 



Moment Order P* Theorem # 4 Bounds EMM Bounds 

6 .934 < Egr < 1.170 .934 < Egr < 1.150 

7 1.021 < Egr < 1.168 1.028 < Egr < 1.153 

8 1.027 < Egr < 1.080 1.028 < Egr < 1.067 

9 1.050 < Egr < 1.068 1.059 < Egr < 1.067 

10 1.055 < Egr < 1.063 1.059 < Egr < 1.062 

11 1.055 < Egr < 1.062 1.059 < Egr < 1.061 

12 1.0602 < E,,,. < 1.0613 1.0602 < E,,,. < 1.0610 



P* : {f^2p\0 <p<P} 
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V. CONCLUSION 



We have outlined a theoretical procedure for transforming Barta's configuration space 
theorem into a moment problem equivalent. The advantages of the latter are that it leads to 
a (quasi)-convexity/concavity reformulation that avoids multi-extrema difficulties associated 
with the configuration space formulation. In addition, by so doing, we solve the problem of 
defining a procedure for improving Barta's bounds, once an initial trial configuration is used. 
This was an outstanding, theoretical problem, within the configuration space formulation. 
We show that the Eigenvalue Moment Method (EMM), is an integral part of this procedure, 
and allows us to prove Theorem # 4. In turn, the results presented here prove that the EMM 
feasibility energy values correspond to a continuous set (an interval) since it is bounded by 
the suprcmum and infimum of the extremal eigenvalues associated with the underlying 
Generalized Eigenvalue Problem. 
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